Model for Dynamic Self- Assembled Magnetic Surface Structures 
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We propose a first-principles model for self-assembled magnetic surface structures on the water-air 
interface reported in earlier experiments [l!,!3- The model is based on the Navier-Stokes equation for 
liquids in shallow water approximation coupled to Newton equations for interacting magnetic parti- 
cles suspended on the water-air interface. The model reproduces most of the observed phenomenol- 
ogy, including spontaneous formation of magnetic snake-like structures, generation of large-scale 
vortex flows, complex ferromagnetic-antiferromagnetic ordering of the snake, and self-propulsion of 
bead-snake hybrids. The model provides valuable insights into self-organization phenomena in a 
broad range of non-equihbrium magnetic and electrostatic systems with competing interactions. 

PACS numbers: 87.19.ru; 43.35.Lf; 75.50.Tt; 47.63.Gd 



Understanding the fundamental principles governing 
dynamic self-assembly in non-equilibrium systems con- 
tinues attracting enormous attention in physics and en- 
gineering communities [IHl]- The interest is stimulated 
by the need for creating smart, dynamic materials capa- 
ble of self-assembly, adaptation to environments, and for 
the design of artificial structures capable of performing 
useful tasks on the microscale [ol fl^ , including targeted 
cargo delivery [11] or stirring in microfiuidic devices [l^ . 

In a series of vi^orks [l|, Q we reported experimental 
studies of self-assembled dynamic magnetic microstruc- 
tures (magnetic snakes) formed on the water-air inter- 
face from a dispersion of magnetic particles energized by 
an alternating (ac) magnetic field applied perpendicular 
to the interface. The snakes form due to coupling be- 
tween the fiuid's surface deformations and the collective 
response of particles to an ac magnetic field. These spec- 
tacular dynamic structures have complex magnetic order- 
ing the snakes' segments are formed by ferromag- 
netically aligned chains of microparticles; however, the 
segments are always anti-ferromagnetically ordered, see 
Fig. [TJ The snakes generate two pairs of large-scale vor- 
tices located at the tails. Under certain conditions snakes 
spontaneously break the symmetry of the vortex pairs 
and turn into self-propelled entities 0, • Some aspect 
of the snakes' behavior were reproduced by phenomeno- 
logical models based on the Ginzburg-Landau type equa- 
tion for surface waves coupled to a large-scale flow 
Nevertheless, the fundamental microscopic mechanisms 
leading to the formation of snakes and the relation be- 
tween properties of the snakes and microscopic properties 
of the system remain unclear. 

In this Letter we report on a mathematical model 
capturing entire self-assembly dynamics in a system of 
magnetic microparticles on the water-air interface. The 
model is formulated in terms of the Navier-Stokes equa- 
tion for fluids in a shallow water approximation coupled 
to interacting magnetic particles described by the corre- 
sponding Newton equations. The particles-fluid coupling 




FIG. 1. (Color online): Self-assembled magnetic snake gener- 
ated by an ac field (/ = 50 Hz, Ho = 100 Oe), composed of 
90/im Nickel spheres: individual particles form self-assembled 
ferromagnetic chains which in turn align in segments of the 
snake. Arrows show the orientation of the total magnetic 
moments in the segments. Inset: Pair of counter-rotating 
vortices at the tail of a magnetic snake. 



is two-fold: an external ac magnetic field causes oscil- 
lations of the particles and deformations of the air-fiuid 
interface. In turn, hydrodynamic flows advcct/reorient 
particles and thus mediate magnetic dipole interactions. 
Our numerical studies faithfully reproduce observed phe- 
nomenology: formation of dynamic snake-like arrange- 
ments of magnetic particles from initially disordered 
configurations, complex ferromagnetic-antiferromagnetic 
ordering of the snake, quadrupole structure of self- 
generated vortex flows, and even self-propulsion of a 
bead-snake hybrid. The computational algorithm is im- 
plemented for graphics processing units (GPUs). Insight 
and techniques developed in the course of this work can 
be applied to a variety of non-equilibrium interfacial non- 
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FIG. 2. (Color online): (a), (b): Surface flows generated 
by dimers/short chains due to external torques exerted by a 
vertical ac magnetic field, (c): coupling between chains and 
surface deformations. White arrows indicate the magnetic 
moment orientations in segments. Ht is a tangential compo- 
nent of the external ac field enforcing a specific orientation of 
the chain's magnetic moment. 



equilibrium systems with competing interactions. 

A precise description of the motion of particles sus- 
pended on the air-water interface is computationally pro- 
hibitive. In particular, accurate modeling of the three- 
dimensional finite Reynolds number (Re) Navier-Stokes 
equation (snakes have Re ~ 100), realistic account for 
particle-fluid interaction and description of solid-fluid 
contact lines required enormous computational power. 
Here we propose a simplified yet nontrivial model al- 
lowing for the investigation of the entire process of the 
snakes' self-assembly. The model identifies key physical 
mechanisms and ingredients. We made the following ma- 
jor approximations: (i) we considered the Navier-Stokes 
equation in the so-called shallow water limit (see, e.g. 
15|). This approximation is valid when the characteris- 
tic length of the snake is large compared to the thickness 
of the fluid layer /iq. We verifled independently in our 
experiments that the snakes exist on thin layers of liq- 
uid, (ii) We simplified the magnetic interaction between 
particles by restricting orientation of their magnetic mo- 
ments along the surface of the fluid. This implies that the 
"elementary" object in our description is a microchain or 
a magnetic dimer formed by two spherical ferromagnetic 
particles (the magnetic moment of the chain is directed 
along the chain) . (iii) We simplified the particle-fluid in- 
teraction by assuming that the primary effect of the ver- 
tical ac magnetic field is a rocking motion of the particles 
leading to the generation of surface flows, see Fig. [2^,b. 
We took into account advection of particles by surface 
flows and rotation of their orientation due to vorticity 
and shear. However, we neglected hydrodynamic flows 
due to in-plane drift and rotation of the particles: for 
the conditions of our experiments, the magnitude of these 
motions is substantially smaller than those produced by 



the external magnetic field. 

The evolution of the fluid surface h and the depth- 
averaged hydrodynamic velocity v is described by 

dth + V(/iv) = 0, (1) 
Qtv + (vV)v = 77(V^v - Cv) - V/i + crVA/i 

+Hq sin(a;i) ^ s(r - rj)pj (2) 

3 

where rj is the kinematic viscosity; a the surface ten- 
sion; the term — describes friction with the bottom of 
a container (we used ^ = 3 for laminar flow). Eq. ^ 
is linearized near the equilibrium height h^. The vari- 
ables are scaled as follows: coordinates r — )• r/ft-o, time 



t t^/ho/g, velocities v — > v/^/gho, where g is the grav- 
itational acceleration. The last term in Eq. ^ describes 
the forcing on the fluid induced by motion of particles in 
the applied ac magnetic field with the magnitude Hq and 
frequency lu. Here, p — (cos((/)), sin(0)) is a unit orien- 
tation vector parallel to the particle's magnetic moment 
and s(r) is a function describing the shape of the particle. 
In most of our simulations we used: s = Ag exp(— t^/sq) 
for r < ro and s — for r > rg, where sq is the parameter 
related to the particle size, rg is a cutoff radius, typically 
5-7 particles sizes, and the parameter Aq (set later to 
Aq = 1) characterizes the strength of the magnetic forc- 
ing. Our simulations show that this particular choice of 
the function s provides a good coupling between particles 
and the continuum hydrodynamic field. 

Particle positions and orientations , pj are governed 
by the following equations [3| 



mpi-j + ^t'i^j = Fj + fit^ — /3V/i 



(3) 



firepj • E • (I —PjPj) X Pj + kHq sin(wt)V/i x pj (4) 

where rrip, Ip, fj,t, and fir are the particle's mass, moment 
of inertia, translational and rotational viscous drag coeffi- 
cients, respectively. Fj and Tj are forces and torques due 
to magnetic dipole-dipole interaction and steric repulsion 
between the particles (torques have only one non-zero 
component), — dvy/dx — dvx/dy is the vorticity and 
Eki — {dvk/dxi + dvi/ dxk) /2 is the rate of the strain 
tensor of the hydrodynamic flow, and I is the identity 
matrix. While we kept - for the sake of completeness 
- the accelerations mprj , Ip(j)j in the equations of mo- 
tion, for typical experimental conditions these terms are 
irrelevant and can be omitted. Eqs. ([3])&(|31) have the 
following meaning: in addition to magnetic forces and 
torques Fj,Tj, the particles are subject to advection by 
the hydrodynamic flow (~ v), sliding down the gradient 
of the surface due to gravity {/3Wh) and rotation due to 
the flow's vorticity (~ fi). Moreover, anisotropic objects, 
such as microchains formed by magnetic particles, will 
tend to align along the axis of the rate of the strain ten- 
sor Eki with a certain prefactor e which in turn depends 
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FIG. 3. (Color online): (a)-(c) Formation of a snake from 225 initially randomly distributed particles. The background color 
represents the height of the fluid surface [h), black arrows depict the orientation of the magnetic moments of the particles, and 
particles are shown as grey circles. Only a part of the entire integration domain is shown, (d) Flows generated by the snake in 
the entire domain, (e) Antiferromagnetic order between the snake's segments; chains of particles are ferrromagnetically ordered 
in each segment. The parameters of the simulations are: amplitude Ho — 0.56, frequency w = 1, layers thickness ho — 1, 
domain area 160 x 160, magnetic moment fio ~ 0.26, viscosity rj = 0.01, and particle diameter a = 0.8, and k — 3 (Movies 
and other parameters are in [l3l)- These parameters are estimated for the following experimental values: layer depth 1 mm, 
particles radius 0.4 mm, frequency of the ac field f ~ 18 Hz (a; = 2nf), and for the viscosity and surface tension of water. 



on the shape of the particle (for spheres e = 0). The last 
term in Eq. ^ deserves a special attention: it describes 
the magnetic alignment of particles along the direction 
of projection of the external ac magnetic field parallel to 
the surface of the fluid, see Fig. [2}:. The in-plane com- 
ponent of the field, oscillating with the same frequency 
Lu appears due to the deformation of the surface, and is 
proportional to V/i. This in-plane field promotes anti- 
ferromagnetic ordering of the neighboring segments and 
is crucial for the formation of snakes. 

The interaction betvifeen particles j, k is described by 
the Hamiltonian H — H'^ + V.'^ , where 'H'' is due to the 
magnetic dipole-dipole interaction. 



[3(Pjejfe)(pfeejfc) - pjpfc] 



(5) 



rjk is distance between particles, e^fe = rjk/rjk, /io 
magnetic moment of individual particle, and H'' = 
fiKa/rjk)^"^ /I6na^ models a sufficiently rigid short-range 
hard-core repulsion between the particles of diameter 
a Correspondingly, forces are evaluated as Fj = 

—dTi/drj and torques as Tj = —pj x dH/dpj. 



Equations ([D)-(|ll) were solved in a periodic x,y do- 
main by a quasispectral method. We used a domain 
area of 160^ in dimensionless units (the length is nor- 
malized by the layer height ho), on a grid with 1024^ 
points, and up to 225 magnetic particles. The algorithm 
was implemented for massive parallel CPUs and run on 
a NVIDIA GTX285 GPU with a peak performance of 1 
TFlop. Typically, a speed up more than 100 times was 
achieved compared to a fast Intel i7 CPU, for details on 
the implementation of the algorithm see [17| . 

Figures[3l^a)-(c) illustrate the formation of a snake from 
initially random distributed particles. The result quali- 
tatively reproduces experimental observations: magnetic 
particles assemble into short chains, then chains form 
segments consisting of several parallel chains, and finally 
segments form snakes-like linear objects. Like in exper- 
iments, the segments are ordered anti-ferromagnetically 
Fig. [31(e). We also plotted the mass flux vector /iv (unlike 
in an experiment, Eq. ([3|) does not provide information 
on the surface velocity). Remarkably, we observed four 
large vortices at the tails of the snake, see Fig. [3] (d). 

We investigated the effect of the amplitude Hq and 
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FIG. 4. Snake's width W vs magnetic field ampfitude Ho in 
simulations (main plot) and in experiment (inset). 
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frequency uj of the external magnetic field on the snake 
structure. Good qualitative agreement was obtained: like 
in experiments (the inset of Fig. 21), we observed a re- 
duction of the snake's width with the increase of Hq, 
see Fig. |4l An increase of the amplitude Hq above a 
certain threshold resulted in the breakdown of a single 
snake and formation of multiple snakes with some aver- 
age length which further decreased with the increase of 
Hq, also in agreement with experiments. An increase of 
the frequency oj resulted in a decrease of the segment's 
size in accordance with the dispersion relation. We also 
studied the effects of variation of the particle's magnetic 
moment hq. We noticed that the snakes are formed only 
in a certain range of /zq. For small fiQ no formation of 
chains was observed, and for too large /io the magnetic 
forces overwhelmed hydrodynamic interactions, and the 
dominant structures were long chains and closed rings Q . 

Our experiment [2] revealed swimmers formed by a 
snake with a non-magnetic bead. The bead at one of 
the snake's tails breaks the balance between vortex fiows 
and turns the snake into a self-propelled entity. To verify 
that, we modified one of the particles in our simulations: 
the last term Hq sin(wt) in Eq. ([2|) for the last particle was 
replaced by —r]C,iv{rj) (i.e. an increased friction around 
this particle with a certain friction coefficient Ci ) ; also the 
size of this particle was increased by a factor of 6 com- 
pared to magnetic particles. The effect of self-propulsion 
of a snake with a bead was successfully captured by our 
model, see Fig. [5l Remarkably, a swimming snake-bead 
hybrid was formed from random initial conditions: the 
non-magnetic bead is expelled to the periphery, spon- 
taneously attaches to one of the tails of the snake and 
forms a swimmer (note asymmetry between the fore and 
aft vortex pairs). 

In conclusion, we developed a microscopic model for 
self- assembled dynamic magnetic structures at the water- 
air interface. We identified the minimal ingredients nec- 
essary for the description of a highly nontrivial process of 
dynamic self-assembly. Our work reveals that the snakes 



FIG. 5. (Color online): Snake-bead hybrid; the non- magnetic 
bead is shown as large black circle, the snake consists of 64 
particles and Ho = 0.65 (other parameters as in Fig. |4]), see 
also Movies 2,3 in [13]. The time difference between frames 
is 85250. (a) and (b) show the same area of the simulated 
system. 

are formed due to a subtle balance between magnetic 
and hydrodynamic forces. The concepts can be applied 
to a wide range of interfacial particle systems driven far 
from equilibrium by external forces, both on the micro 
and nano scales. The research was supported by the U.S. 
Department of Energy, Office of Basic Energy Sciences, 
Division of Materials Science and Engineering, under the 
Contract No. DE AC02-06CH11357. 
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